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We present a mesoscopic hydrodynamic description of the dynamics of colloidal suspen- 
sions. We consider the system as a gas of Brownian particles suspended in a Newtonian 
heat bath subjected to stationary non-equilibrium conditions imposed by a velocity field. 
Using results already obtained in previous studies in the field by means of a generalized 
Fokker-Planck equation, we obtain a set of coupled differential equations for the local dif- 
fusion current and the evolution of the total stress tensor. We find that the dynamic shear 
viscosity of the system contains contributions arising from the finite size of the particles. 

PACS numbers: 

I. INTRODUCTION 



Due to the theoretical and industrial importance of colloidal suspensions, the rheological properties 
of these systems have been the subject of intense investigations by using a variety of experimental 
and theoretical methods 

mm 

. Among others, the experimental techniques used to study 
these systems are neutron and light scattering, and ultrasonic absorption 1^, Q, 0]. Frorr^t 
theoretical point of view, non-equilibrium statistical mechanics and thermod^ 



10|, 



11 




1^ . mode-coupling formalism 13] and computer simulations 
techniques in which particular interest has been put on the behavior of the frequency-dependent 
shear viscosity of semidilute suspensions of hard spheres in different situations involving direct or 



hydrodynamic interactions 



In this article, we present a mesoscopic hydrodynamic description of the dynamics of colloidal 
suspensions subjected to stationary non-equilibrium conditions imposed by a velocity field. In 
particular, we analyze the behavior of the dynamic viscosity due to its importance in reflecting the 
microstructure details into the relaxation of the system under the conditions imposed by the applied 
flow. At mesoscopic level, the contributions of the suspended phase to the dynamics of the whole 
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system can be described by using a Fokker-Planck equation |lCll.ll9l.l20ll2ll|. This equation has been 
previously obtained on the grounds of irreversible thermodynamics in the dilute regime 3, 
However, it is important to stress that similar equations have been derived in the context of different 



approximations ranging from the kinetic theory of gases to the projector operator formalism 
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22] • From this Fokker-Planck equation we calculate a set of hydrodynamic equations for the 



moments of the distribution function. Then, by using the fluctuating hydrodynamics approach we 
obtain a frequency dependent correction to the dynamic viscosity of the system that accounts for 
memory effects |23l . |24| . 

The paper is organized as follows. In Section II we summarize the derivation of the Fokker- 
Planck equation and carry out the hydrodynamic description of a gas of Brownian particles. In 
Section III we derive the dynamic viscosity of the system and present numerical results. Finally, 
Section IV is devoted to conclusions. 



II. KINETIC AND HYDRODYNAMIC DESCRIPTION OF A GAS OF BROWNIAN 

PARTICLES 

We will analyze the dynamics of a system consisting in a "gas" of Brownian particles of mass m 
suspended in a newtonian heat bath subjected to stationary non-equilibrium conditions imposed 
by a stationary flow Vo{f). 

The usual starting point to describe the dynamics of this system is the A^— particle distribution 
function and its corresponding multivariate Fokker-Planck equation. In this equation, microscopic 
expressions for the hydrodynamic and direct interactions must be introduced. In the semidilute 
regime, for example, one must introduce the Oseen tensor and the two-particle interaction potential. 
At higher concentrations, the subsequent approximations may, in general, consider the Rotne- 
Prager tensor and the three-particle (or higher) interaction potential. It is clear that this approach 



constitutes a A^— body problem hard to solve by analytical methods |25l . [2 

Here, we will follow an alternative approach to the problem in terms of a mean field approxima- 
tion in which the difficulties of the A^— body problem are avoided. This approximation, based on 
the one-particle distribution function, allows one to make some analitycal progress when calculating 
the hydrodynamic equations, despite that the explicit dependence of the transport coefficients on 
frequency and wave vector and the form of the direct interactions are a priori unknown. However, 
by assuming in first approximation that the transport coefficients appearing in the hydrodynamic 
equations are constants, we may use of the well established formalism of fluctuating hydrodynamics 
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to obtain the explicit expressions for the corrections to the transport coefficients as a functions of 



the frequency and the wave vector 2^ 
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According to this, at mesoscopic level the dynamics of the dilute Brownian gas can be described 
in terms of a single-particle local probability distribution P{f, u, t) depending on time t and the 
instantaneous position and velocity of the Brownian particle r and u. This distribution function 
obeys the continuity equation 

|p + V.(SP) = -|=.(PV'.), (1) 

where V^g is a streaming velocity in n— space. The explicit form of V^g can be obtained from the 
entropy production of the system, after assuming valid the hypothesis of local equilibrium in phase- 
space and using the rules of mesoscopic non-equilibrium thermodvnamics llfll . l28i| . Then, entropy 
production can be calculated from the nonequilibrium entropy functional (271] 

ps{f, t) = -ks / Piu, r, t) In —du + psie{f), (2) 

where ks is Boltzmann's constant and ps{f,t) and psi(,{r) are the nonequilibrium entropy per 
unit volume and the local equilibrium entropy per unit volume, respectively. Moreover, we have 
introduced the fugacities z = P{f,u,t)A{p) and zie = exp |[//ze ~ ^ ~ vo)'^]m/kBT^ with pie the 
chemical potential at a local equilibrium and p{f, t) = m J P{u, r, t)du the mass density field of the 
Brownian gas. Finally, A{p), related with the excess of osmotic pressure of the system, represents 
the activity coefficient that introduces the interactions among Brownian particles [^l- 

After calculating the entropy production using @ , identifying forces and currents and using the 
linear-law scheme of non-equilibrium thermodynamics 29] , one finds that the distribution function 
P{r, u, t) obeys the Fokker-Planck equation 



+ V • (uP) - V 



^^^lnA(p) 



m 



dP _ d 
du du 



m ou 



(3) 



where ksT is the thermal energy (assumed constant). In the general case, the friction tensor 
l3{r,t) must contain the average effect of hydrodynamic interactions. It is related to the Onsager 
coefficient i{r,t) through (5 = ^ + e ■ ViTo) in which the Onsager coefficient e{f,t) accounts for the 
effects of the force exerted by a fiuid on a particle of finite size according to the Faxen theorem 
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It is important to notice that vq is the stationary solution of the Navier-Stokes equation 
PqVqVvq = —V • P*^ for the heat bath, with the appropriate boundary conditions on the surface of 
the colloidal particle and pQ is the constant density of the bath. 
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Equations similar to (131) have been derived in the context of different approximations, see for 
example, Refs. Q, H, H 



221 



From equation Q, one may calculate a set of evolution 



equations for the moments of the probability distribution, giving the "hydrodynamic" description 
of the system 2, . After taking into account the usual definitions for the hydro dynamic-like 
fields: mass density p{r,t), velocity v{f,t) and pressure tensor P'^(r, t) (see reference |ld|), the first 



equation of this set corresponds to the conservation of mass 



where the diffusion current JD(r, t) = pv of the Brownian gas satisfies the evolution equation 



27 



2t\ 



Jd+P 



— 1 



DJp 
Dt 



pvo 



m 



dlnp 



V-P^(f,t) 



(5) 



Here D/Dt is the material time derivative, /c^ is Boltzman's constant, and T is the temperature 
of the heat bath. The matrix of relaxation times of the diffusion current is: r i ~ /3 . Finally, Eq. 
© contains the pressure tensor F\f, t) of the suspended phase whose evolution equation is [lO| 



m 



-P 



(3 - e- Vvo 



1 D 

2m^ 



(6) 



where the symbol s stands for the symmetric part of a tensor and T2 is a matrix of relaxation times 
given by 

-1 



(7) 



where we have neglected inertial terms proportional to V • u, assumption valid for small density 
fluctuations Ap/p <^ 1. 

The explicit expression for the pressure tensor can be obtained by solving the set of differential 
equations ©• However, in order to make analytical progress, let us assume the approximation 
1 — /? • \/v I with 1 the unit tensor, and which allows to write Eq. (jSj) in the form 



T2 



— 1 

~ (3 



ni 



pi -\Dop- 



i-i\t 



/3 



^1 

Dt 



(8) 



where the f represents the transpose, we have defined the diffusion coefficient Dq = ^p/? and 
assumed that the leading term in the matrix of relaxation times is (3 . Now, the assumption of 



small density fluctuations implies that the diagonal terms of P'^ relax faster than the non-diagonal 
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ones. Therefore, Eq. © can be separated in the form P*^ 
stress tensor. This separation leads to the equations 



+ L, with L(f, t) the traceless 



„id 



P 



m 



-P 



and 



2L + 



i\ t 



2\Dop- 



After substitution of equations © into © we obtain 



Jd = pvo - D{p) -Vp-p • (V • L 

where we have defined the coUective diffusion coefficient 

din A{p) 



—1 



DJp 
Dt ' 



D{p) = Do 



in accordance with Ref. 



1 + 



dlnp 



Dq—^ . 

op 



(9) 
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29|. Notice that in the hmit of very dilute suspensions, the activity 
coefficient vanishes and then Eq. ()1U() recovers the expected form. 

The evolution equations governing the behavior of the hydrodynamic fields describing the whole 
suspension as an effective medium can be obtained from Eqs. Q and pOj) by adding the local 
diffusion current J = Jo — pvQ and the traceless part of the stress tensor of the heat bath given by 
P'' = — 2r/o (Vuo)'^ with rjQ the Newtonian viscosity of the bath 36 1. 

Thus, multiplying Eq. Q by pop~^, adding P*', neglecting the convective term in the total 
time derivative and using the definition of J, one arrives at the evolution equation for the total 
traceless stress tensor of the suspension Q'^ = P'^ + Ipop^^Ij 



dt 



Q 



1 ° 









= -2 


DoPoPb^ • (VJ) 


- 2 


Vfs ■ (Vwo) 



(12) 



where pB = rnn, being n = N/V the number density of colloidal particles, and we have defined 
the effective viscosity coefficient fj of the whole suspension as 



Vfs=m 



(13) 



Equation (|13|) contains corrections due to the finite size of the particles through the coefficient 
e ~ (m/G/csT) • /3 with a the radius of the particle. The expression for e has been obtained by 



28i |. In deriving Eqs. (|12() and ((T^ we have assumed that 



taking into account the Faxen theorem 
terms of the order DqJp^^ ■ Vp may be neglected and that the coefficient of the first term at the 



right-hand side of (|T2|) only involves the average constant value of the density of the suspended 
phase. This is consistent with the approximation Ap/p <^ 1 used to derive Eq. ©. 

Now, by substituting Jo = J + pvo into Eq. (fTU]) and using the Navier-Stokes equation, one 
can derive the evolution equation for the local diffusion current J 

DJ 



qs- 



(14) 
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where we have introduced the effective diffusion tensor 



D,ff = D{p)-^^ \Do-\^[{l + f)-Vvo]\p'\VJfY (15) 

and Q^g is the quasi-stationary form of the stress tensor of the suspension whose exphcit form 
follows from Ea. (ll2l) for times t ^ (/3~^)/i- The dependence of Dpff on the stress tensor of 

inn 

the heat bath has been previously reported in Refs. |1(]| . 12(1 . 122l |. Now, in accordance with the 
assumption of low velocity gradients, the term (j ■ Vwo can be neglected in Eqs. (|14|) and (|15|) . 
and then we obtain the simplified equation 

—1 DJ - 



(16) 



where we have assumed that Q'^ can be introduced instead of Q^^, and that the coefficient of the 
last term only involves the average constant value of the density of the suspended phase pB- 

Equations similar to 1)12(1 and have been discussed by several authors from different macro- 
scopic and microscopic points of view in the context of polymer solutions 0,0, These 
equations incorporate viscoelastic effects into the description since they contain a term associated 
with the relaxation of the involved quantities (Q*^, J), and thus lead to more general constitutive 
relations than those for Newtonian fluids. Equations (jl2(l and ()l(i() include tensor coefficients being 



then suitable to describe anisotropic systems 



35l |. In the following section, we will use equations 



(|12|) and (|16|) to analyze the rheological properties of the system when fluctuations of the suspended 
phase are considered. 



III. THE DYNAMIC VISCOSITY 



The dynamic viscosity coefficient of the suspension as a function of the frequency uj, can be 
calculated by considering the linearized fluctuating equations for the hydrodynamic fields which 
follow from equations Q), (|T^ and (fT?T|) . With this purpose, we will assume that p = ps + Sp, 
vo = Vs + Su and Q*^ = Qq + SQ^, with 6p, 6u and 6Q^ the deviations with respect to the average 
values pbi Vs and Qq. 

Moreover, in this section we will consider the case of a dilute suspension of particles. This 
assumption leads to a simplification of the description in which both hydrodynamic and direct 
interactions can be neglected. As a consequence, the set of Onsager and transport coefficients be- 
come scalar and, in the simplest approximation, constant: /? = e*= el, = ^1 and fjj^ = rffgl, 
and the activity coefficient A(p) = 1, thus simplifying the expression for the diffusion coefficient to 
D{p) = Do. 
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Then, at first order in the deviations, from Eqs. ()12|) and (|16j) one obtains the fohowing 
set of fluctuating hydrodynamic equations 

d{6p) 



dt 



-V • ((5J) -Vpo-6u-V {6p) ■ Vs, (17) 



6J + P~^§i^J = -DoV6p + p-'pBPo'y ■ SCf, (18) 



<5Q° + /3-^-^5Q° = -2Z)oPBP^' (V5J)° - 2??/, {V6uf . (19) 

To obtain the corrections to the dynamic viscosity of the suspension, Eq. H19|) can be written in 
a more suitable form by taking into account the separation vq = Vs + Su, which can be interpreted 
as the sum of the unperturbed imposed velocity field Vs , in particular, a shear flow, and 5u the 
corresponding correction to the velocity field arising from the presence of the particle, which is 
proportional to Vvs 30]. Following Ref. |2fl], Eq. ^ can be written in the form 

5Q° + /?^';^'^Q° = -2DopBPo^ (V6jy - 5<P7]fs (Vvs )° . (20) 

where (p represents the volume fraction occupied by the Brownian particles in the whole system. 

Now, by taking the Fourier transform of Eqs. (|17j) . (|18jl and (|20p and solving for 5Q,^ in terms of 
ik6u, and taking into account that the frequency-dependent correction to the viscosity Atj (Jc, lo^ 
of the suspension is defined according to the relation 

(5Q° = -i[Ari (fc,a;) k5u]^, (21) 

we find that 

Ari (k,uj)c^ -<P "^f' n . 1.2 , (22) 



where, for simplicity, we have considered the case in which k = (/!;,0,0). Notice that the frequency 



dependence of the c 
fraction expansion 



namic viscosity for semidiluted suspensions follows a second order continued- 



ynam 

i If 



we now take into account the relation Q'^ = Qq + 5Q^, we obtain that 



the total viscosity of the suspension rjfg + Arj {k, to) can be written as 

5 . 1 



r]{k,uj) ~ r]fs 



1 + 



(23) 
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It is interesting to notice that by expanding Eq. (|23() up to first order in to and second order in k 
and taking into account Eq.(^Sl)) we obtain 



5 , im ,oDo 



1 + Sc ^ (1 + e)l , (24) 



where Sc = f]o/{poDo) is the Schmidt number for the Brownian gas diffusing in the solvent. Eq. 
(|24|1 is an expression similar to that reported in Ref. 3 a ]. Nonetheless, Eq. (|24() includes corrections 
due to the finite size of the particles which are 



usme 



To analyze the frequency behavior of Eq. (|23|) . it can be recasted in a more convenient form by 
the dispersion relation k'^ = — z^^^, which can be obtained from the Navier-Stokes equations 



28 



by 



30|. Substituting the last expression into Eq. (pS)) and using dimensionless variables defined 



^^^^ k^^^k' and (25) 
we obtain the following expression for the dynamic viscosity 

+ i ^ , (26) 

bc[l—iui) 

which only depends on the dimensionless frequency a;, the volume fraction and the Schmidt 
number Sc. The term containing Sc represents a contribution arising from the coupling between 
diffusion and viscosity. Note that our definition of the Schmidt number is related to that used 
in Ref. [s^l by means of the relation Sc = ScDj/Dq, being Dj the fluid self-diffusion coefficient. 
For a gas Sc ^ 1 (this value is also used in numerical simulations of colloidal systems as, for 

n 

example, the one of Ref. 15]), whereas for a liquid Sc^ 1 and since for a colloidal suspension in 
general Df ^ Dq then Sc ^ 1. Therefore, the terms containing Sc ^ in Eas. (|13() and (|26() can_be 



neglected and the expression for the dynamic viscosity reduces to a form similar to that of Ref. 
the differences being the finite size contributions contained in rjfg, the explicit dependence on the 
volume fraction and that the characteristic relaxation time is given by instead of td = a? / Dq. 

In Fig. 1 we show the real rj' (cj) and imaginary rj" (cj) parts of the complex dynamic viscosity 
r/ (tj) as a function of logcj for different values of the volume fraction and considering Sc ^ \. 
Notice that ri{u)) has the usual behavior of a relaxation processes with only one relaxation time 
^. From Ea. ()26() it can be seen that the ratio [ry' (lj) — 1] / [??o — 1] , with ?7g = i]' (cj = 0), will 
depend only on the dimensionless parameter cj. The imaginary term can be reduced in the same 
manner rj" [to) / [r/g — 1]. In Fig. 2 we plot the reduced viscosity as a function of to. In addition to 
the usual case of very large Schmidt numbers (solid line), we have also explored what would be the 
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consequences when the Schmidt number is of the order of unity. For Sc= 1 (dashed hne) the real 
part of the viscosity presents a hump near uj = 1. This characteristic appears for values of Sc in 
the range [1,2 + \/3] whereas for Sc>2 + -v/3 the hump is not present (for Sc> 1, the imaginary 
part of rj (uj) takes only positive values). At intermediate frequencies rj' {uj) decays following the 
law u}~'^, independently of the value of Sc. This is a typical behavior for processes with only one 

imaginary part of the viscosity has a bell-shaped form 
for 5*0 S> 1, whereas for Sc = 1 it becomes asymmetric and more stretched. It is important to 
notice that for intermediate frequencies, rj" {u) decays as also a usual behavior for processes 
with only one relaxation time. The low frequency region of the dynamic viscosity curves deserves 
also a detailed analysis. In this frequency region, rj" {uj) grows as tD, independently of the value of 
the Schmidt number, except for Sc = 1 where the curve grows as a;^. 

Finally, it is important to point out that it has been experimentally observed that for interme- 
diate frequencies both, the real and imaginary parts of the viscosity decay as ijJ~^^'^ Such a 
frequency behavior can be modeled by assuming a distribution of relaxation times Tp. In this case, 
the viscosity can be written as 

P=l ^ Sc{l-lOJTp) 

where N is the number of relaxation processes involved. Gp is a constant relaxation strength of 
the relaxation process with relaxation time Tp, given by a distribution of the form 

Tp = Tip-"", a = 2 (28) 

with Ti the longest relaxation time 0]. For a non- interacting hard-sphere system this would 
correspond to a distribution on the size of the particles Up = aip", or equivalently, Scp = Scip". 
With this size distribution our results for the viscosity lead to the same power law decay on 
frequency a;~^/^. In Fig. 3 we show this behavior for two different values of the Schmidt number 
corresponding to the longest relaxation time ti, Sci = 1, and Sci ^ 1. We notice that the 
intermediate frequency power-law dependency on frequency is not altered by this choice but at 
low frequencies differences shows up for fj" (uj). A different power law for the size distribution will 
modify the behavior of the viscosity. 



IV. CONCLUSIONS 



We have proposed a model for the dynamics of colloidal suspensions based on a mesoscopic 
hydrodynamic description derived by applying the rules of mesoscopic non-equilibrium thermody- 
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namics. We have derived the viscosity of a monodisperse hard-sphere system that includes effects 
due to the finite size of the colloidal particles. It has been shown that the normalized dynamic 
viscosity depends only on three dimensionless parameters: the Schmidt number Sc, the volume 
fraction 0, and the normalized frequency uj. Also, it presents power-law exponents for different re- 
gions of the frequency space. Our results show that the mesoscopic approach leads to results which 
are consistent with previous theoretical descriptions derived by using the generalized hydrodynam- 
ics theory. We expect that the present model may be extended to a wider class of suspensions 
including particles with internal structure. 
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FIG. 1: (a) Real and (b) imaginary parts of the normalized dynamic viscosity as a function of a; for Sc^\. 
Volume fractions as indicated in the figures. 
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FIG. 2: (a) Real and (b) imaginary parts of the reduced dynamic viscosity as a function of w. The solid 
line corresponds \.o Sc^ \ and the dashed line corresponds to Sc= 1. 
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FIG. 3: (a) Real and (b) imaginary parts of the reduced dynamic viscosity as a function of ui, for a 
distribution of relaxing processes. The solid line corresponds to Sci 3> 1 and the dashed line corresponds 
to Scx = 1. The number of relaxing processes considered was N = 1000. 



